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We have performed realistic shell-model calculations for nuclei around doubly magic i™Sn and 
^''^Sn using an effective interaction derived from the Bonn A nucleon-nucleon potential. The results 
are in remarkably good agreement with the experimental data showing the ability of our effective 
interaction to provide an accurate description of nuclear structure properties. 
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I. INTRODUCTION 



■ A fundamental goal in nuclear structure theory is to understand the properties of complex nuclei in terms of the 
(SJ ^ nucleon-nucleon [NN) interaction. Since the pioneering work of Kuo and Brown 0], who in the mid 1960s derived 

an s-d shell effective interaction from the Hamada- Johnston potential Q|, there has been considerable progress in this 
field. On the one hand, the theoretical framework in which the model-space effective interaction Voft can be derived 
from a given {NN) potential has been largely improved (the main aspects of this derivation are reviewed in Ref. Q). 
On the other hand, high-quality NN potentials have been constructed which give a very good description of the NN 
, scattering data. A comprehensive review of modern NN potentials suitable for application in nuclear structure is 

■ given in Ref. 

. These improvements have brought about a great deal of renewed interest in realistic shell-model calculations. In 
OO ■ this context, the main question is to which extent modern realistic interactions can provide a consistent and accurate 
description of nuclear structure phenomena. To try to answer this question, the study of nuclei in the vicinity of closed 
. shell is extremely important. In fact, they provide the best testing ground for the basic ingredients of a shell- model 
' calculation, in particular as regards the matrix elements of the effective NN interaction. 
'Tx i In this paper, we shall present some results of realistic shell-model calculations for nuclei around doubly magic 
^ • loogj^ i32gj^^ rpj^gy ^j^g iv = 50 isotones ^^Cd, ^^Ag, ^^Pd, the TV = 82 isotones ^^^Te and and the doubly 
odd nucleus ^'^^Sb, which has a single proton outside the Z = 50 closed shell and a single neutron hole in the closed 
= 82 shell. These nuclei lie well away from the valley of stability and only recently more experimental information 
has become available which is of great value for a stringent test of our calculated effective interaction. It should be 
noted that the study of ^'^^Sb provides a direct test of the effective interaction matrix elements with isospin T = 0. 
$H ' In earlier works we have performed shell-model calculations for the light Sn isotopes making use of two effective 
5^ , interactions derived from different free NN potentials. They are the Paris |^ and Bonn A |^ potentials. It turned 
out that the latter, which has a weaker tensor force component, leads to the best agreement with experiment for all 
of the nuclei considered For this reason, we have used it in the present work. 

Our presentation is organized as follows. In Sect. 2 we give a brief description of our calculations. In Sect. 3 we 
present our results and compare them with the experimental data. Sect. 4 presents a summary of our conclusions. 



II. OUTLINE OF CALCULATIONS 



As already mentioned in the Introduction, we make use of a two-body effective interaction derived from the meson 
theoretic Bonn A potential. This was obtained using a G-matrix folded-diagram formalism, including renormalizations 
from both core polarization and folded diagrams. A description of the derivation of the effective interaction VcS from 
the nucleon-nucleon potential including references can be found in Ref. [|j . We only outline here the essential of the 
method and point out the main differences between our present and earlier calculations [0-^. 

The effective interaction can be schematically written pO[ | in operator form as 

V,fi^Q-Q' [ Q + Q' [ Q [ Q-Q' [ Q [ Q [ Q ■■■ , (1) 
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where Q and Q' represent the Q-box, composed of irreducible valence-Hnked diagrams, and the integral sign represents 
a generalized folding operation. We take the Q-box to be composed of G-matrix diagrams through second-order in G; 
they are just the seven first- and second-order diagrams considered by Shurpin et al. To calculate the effective 
interaction given by Eq. (1), a first step is to calculate the irreducible Q-box diagrams and their energy derivatives 
in terms of the model-space G-matrix defined by 

G{lo)^V + VQ2 ^—^Q2G{lo). (2) 

Here Q2 is the Pauli exclusion operator for the two interacting nucleons, V represents the NN potential, T denotes the 
two-nucleon kinetic energy, and lu is the so-called starting energy. We employ a matrix inversion method to calculate 
the above G matrix in an essentially exact way . 

For the = 50 isotones, we have considered the doubly closed ^°"Sn as an inert core and treated protons as valence 
holes. This leads to a calculation of the Q-diagrams which is different from the usual one for particles. A detailed 
description of the calculation of our two-hole effective interaction will be given in a forthcoming publication. We 
have chosen the Pauli exclusion operator Q2 in the G-matrix equation (2) as specified [|2| by (ni,n2,n3) = (11, 21, 
45). For the shell-model oscillator parameter Ticu we have used the value 8.5 MeV, as obtained from the expression 
hu; = 45^-1/3 _ 25A-2/3 for A = 100. 

For the = 82 isotones the valence-proton and -neutron orbits outside the ^■^^Sn core are different. Therefore, we 
have chosen (ni, 77.2,713) = (16, 28, 55) for the neutron orbits and (711,712,713) = (11, 21, 55) for the proton orbits. 
The effective interaction has been calculated using an isospin uncoupled representation, where neutrons and protons 
are treated separately. The adopted value of oscillator spacing is hui = 7.88 MeV. It should be mentioned here that in 
our earlier study |^ we made the choice (711,712,713) = (11, 28, 45) for both protons and neutrons. As a consequence, 
the Pauli exclusion operator was not treated in a completely correct way. 

As regards ^^^Sb, we have treated the odd proton and the remaining 31 neutrons as valence particles. This makes 
the T ~ I matrix elements of the effective interaction be just the same as those used in our earlier study of the Sn 
isotopes 0, to which we refer the reader for details. 

For the iV = 82 isotones and for ^'^^Sb our model space includes the five single-particle (s.p.) orbits 0(77/2, 1^5/2, 
2si/2, 1^^3/2, and 0/iii/2, while for the iV — 50 isotones the proton holes are distributed in 0179/2, lpi/2, IP3/2, f^nd 
O5/2 orbits. 

As regards the choice of the s.p. energies, we have proceeded as follows. For the iV = 82 isotones we have taken 
three s.p. spacings from the experimental spectrum of ^^'^Sb |jl^,|I^. In fact, the 37/2, ^5/2, (^3/2 and /111/2 states can 
be associated with the ground state and the 0.962, 2.439 and 2.793 MeV excited levels, respectively. It should be 
noted that the position of the 1^3/2 state corresponds to the new value very recently provided by the high-sensitivity 
7 spectroscopic measurement of Ref. |l^. This is about 400 keV smaller than the previously accepted value (2.708 
MeV) jl^. It turns out, however, that this change has little effect on the spectra of ^^''Te and ^^^I. As for the Si/2 

state, its position has been determined by reproducing the experimental energy of the level at 2.150 MeV in -'^'^^Cs. 
This yields the value es^^^ = 2.8 MeV. 

For the = 50 isotones, the single-hole energies cannot be taken from experiment, since the single-hole valence 
nucleus ^^In has not yet been studied. Therefore we have determined them from an analysis of the low-energy spectra 
of the isotones with A > 92. The adopted values (in MeV) are: eg^^^ — 0-0, Cpi^a = 0.73, Cp^^^ — 2.17, e/^^^ = 3.24. 

Regarding ^'^^Sb, we have assumed the s.p. energies to be the same for neutrons and protons. Our adopted values 
(in MeV) are: e^j.^^ — 0, eg^^^ = 0.20, e^^^^ = 1-72, Cd^^^ — 1.88, and e/j^^^^ = 2.70. As compared to the set of s.p. 
energies used for the Sn isotopes (see Ref. 0), only Csi/a ^^"^ ^<i3/2 have been modified. More precisely, they have 
been both decreased by about 0.5 MeV. It should be noted that the position of these two levels plays a minor role 
in the calculations for the light Sn isotopes while it is very important to satisfactorily reproduce the experimental ij 
state and to place in the right energy range the calculated negative-parity states. 

III. RESULTS 

A. A*' = 50 isotones 

In Fig. 1 we compare the calculated 2+, 4+, 6+ and 8+ yrast states of ^*Cd with those recently identified in the 
study of Ref. [ p^ . We see that our results are in very good agreement with experiment. A measure of the quality of 
the agreement between theory and experiment is given by the rms deviation a uM, whose value is 105 keV. 
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It should be noted that the predicted position of the 5~ state (2.73 MeV) is quite consistent with the experimental 
information available for the two lighter even isotones. In fact, in ®^Pd and ^^Ru a 5~ state has been observed 18 
at 2.648 and 2.625 MeV, respectively. As regards the structure of the excited states, we find that the positive-parity 
states are of practically pure {Trgg/2)~^ character while the 5~ state is dominated by the .9(^/2^1/2 configuration. On 
the contrary, the ground-state has a significantly mixed wave function, the percentage of configurations other than 
(7r5g/2)~^ being about 17%. 

In Figs. 2 and 3 the calculated spectra of ^^Ag and ^^Pd are compared with the experimental ones [|l^,^. We 
see that for both nuclei our calculations produce a spectrum very close to the experimental one. The main points of 
disagreement are the position of the state in ^^Ag and that of the 2+ state in ^^Pd which lie 213 and 281 keV, 
respectively, above the observed ones. The a value is 121 keV and 130 keV for ^^Ag and ^^Pd, respectively. 

Experimental information on electromagnetic transition rates in heavy TV = 50 isotones is very scanty. It is of 
interest to mention here, however, that in the work of Ref. [[| a value of 0.44^°:^° W.u. for the B{E2; 8+ 6+) 
in ^^Cd has been reported. Using the bare proton charge Cp = 1 we obtain 0.67 W.u. By contrast, to reproduce the 
experimental value of the same B{E2) in ^^Pd (0.34 ± 0.05 W.u.) an effective charge of at least 1.7e is needed. To 
clear up this point more experimental data are required. 



B. i^^Te and i^Sj 

The experimental [^,0 and theoretical spectra of the two-proton nucleus ^'^^Te are compared in Fig. 4, where 
all the calculated and experimental levels up to 3.2 MeV excitation energy are reported. We see that while the 
theory reproduces all the observed levels it also predicts the existence of a 3+ and a 0+ state at 2.65 and 2.78 MeV, 
respectively. This prediction is strongly supported by the experimental information available for the two heavier even 
isotones. In fact, in ^^^Xe a 0+ state has been observed at 2.58 MeV while in ^^^Ba both a 0+ and a 3+ state 
have been located at 2.34 and 2.45 MeV , respectively. Above 3.2 MeV excitation energy the comparison between 
theory and experiment is made only for those observed levels which have received a spin-parity assignment. We do 
not include the new levels observed in ||2^ since all of them should be interpreted as neutron particle- hole states. This 
interpretation is confirmed by our calculations. In fact, the only two states, having J'^ = 8"*" and 10^, which can be 
constructed in our model space are both predicted to lie at about 7.3 MeV while the two experimental states with 
these spin-parity assignments have been located at 4.557 and 5.622 MeV, respectively. We see that the calculated 
spectrum reproduces very well the experimental one, the discrepancies between theory and experiment being smaller 
than 50 keV for several states (9 out of 15). The rms deviation a is 106 keV. 

In Fig. 5 we compare the calculated spectrum of ^'^^I with the experimental one [^4[|25| ] up to 4.0 MeV excitation 
energy. As in the case of ^^''Te, we exclude the experimental levels above 4.2 MeV, which originate from core 
excitations. We should note that the spectra of Fig. 5 include all experimental and calculated levels up to 1.5 MeV. 
Above this energy several other levels without assigned spin and parity are reported in [ p5[ ; we compare our calculated 
states only with those observed in . From Fig. 5 we see that the excitation energies are remarkably well reproduced 
for all the reported states, the a value being 58 keV. We have associated the experimental level at 1.010 MeV with 
the theoretical at 0.931 MeV. 

We should now point out that the theoretical results presented here for ^'^'^Te and ^^^I are in a substantially better 
agreement with experiment than those obtained in our earlier study The reason for these improvements can be 
clearly traced to the better treatment of the Pauli exclusion operator Q2 ■ 

In Table I we compare the experimental reduced transition probabilities in ^^''Te with the calculated ones. We 
have used an effective proton charge e^^ — 1.55e. This is consistent with the values adopted by other authors 
1^,^. The theoretical B{E2) values are in very good agreement with experiment. As regards the transitions, 
we find that the _B(£'3; 9]^ 6^) is well reproduced while the B(i?3;9^ is underestimated by a factor of 

about 4. A possible reason for this discrepancy lies in the fact that only a small amount of configuration mixing is 
present in the calculated 6"*" states. In fact, the decay to the 6^ state is dominated by the single-proton transition 
{^11/297/2)9- ^ {97 /2db/2)&+ while that to the Gi^ state by the transition (^11/237/2)9- ^ (37/2)6+1 which is retarded 
owing to spin flip. The theoretical B{E3; 9^ 6'^) value could be brought in agreement with experiment by an 
amount of configuration mixing of about 15%, which would, of course, reduce the B{E3; 9j~ 6^) value. It should 
be noted, however, that the latter would still be within the error bar. 
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C. "^Sb 



The experimental and theoretical spectra of ^'^^Sb are compared in Fig. 6, where all the observed levels are reported. 
In the calculated spectrum all levels up to 1.4 MeV excitation energy are included while in the higher energy region 
only the and 3^ states are reported. It should be noted that the nature of the presently available experimental 
information is quite different for positive- and negative-parity levels. In fact, while the spin-parity assignments to the 
former have been clarified by the study of Ref. this is not the case for the latter. More precisely, the excitation 
energy of the 8~ state is not known (the work of Ref. places it between 150 and 250 keV) and the three other 
observed negative-parity states have not received firm spin assignments. We find that the first excited 8^ state lies 
at 126 keV while the 6j~, 5j~, 7^, 3j~, and 4j~ states are grouped in a very small energy interval (from 210 to 380 
keV). As a consequence, any attempt to establish a one-to-one correspondence between the observed levels and those 
predicted by our calculation could be misleading. It is of interest to note that the above states, which all arise from 
the T^g^ /2^hii/2 configuration, are well separated from the other two members of the multiplet, i.e. the 9~ and 2~ 
states, which are predicted at 1.0 and 1.42 MeV, respectively. A similar behavior is also predicted for the Trdr,/2vh^-^^^ 
multiplet. In fact, the 7^, 6~, 5~, and 4^ states belonging to this configuration lie between 0.82 and 0.97 Mev while 
the highest- and lowest-spin members (8~ and 3~) are at 1.12 and 1.56 MeV, respectively. 

From Fig. 6 we see that the experimental excitation energies of the positive-parity states are remarkably well 
reproduced by the theory, the largest discrepancy being 77 keV for the 5^'" state. The value of the rms deviation <t is 
only 32 keV. 

In Table II we compare the BE{2) values for transitions between states below 1.1 MeV excitation energy with the 
calculated ones. A more complete analysis of the electromagnetic properties of ^^^Sb may be found in js^l- We have 
used an effective proton charge = 1.55e, which is just the same as that adopted for the N = 82 isotones. No 
effective charge has been attributed to the neutron hole. As we see from Table II, the experimental data are affected 
by large errors. In view of this, the agreement between theory and experiment can be considered quite satisfactory. 
In fact, our calculated values lie all but two within the limits set by experiment. 

IV. CONCLUDING REMARKS 

We have presented here some recent results of realistic shell-model calculations for nuclei around doubly magic -'^'^"Sn 
and ^"^^Sn. They have been obtained by employing an effective interaction derived from the Bonn A nucleon-nucleon 
potential. We have shown that the agreement between theory and experiment is very good for all nuclei considered. It 
is to be emphasized that the study of ^'^^Sb provides a test of our T = effective interaction in this mass region. This 
is of special interest since in earlier works using different NN potentials it turned out that not enough attraction was 
provided by the calculated matrix elements of the T = effective interaction, which has a stronger dependence on the 
tensor force strength than the T = 1 interaction (a detailed discussion of this important point including references is 
given in Ref. Q). 

In conclusion, the success achieved by our calculations shows that our effective interaction derived from the Bonn A 
NN potential is able to describe with quantitative accuracy the spectroscopic properties of nuclei near closed shells. 
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TABLE L Calculated and experimental B{E\) values (in W.u.) for ^^''Te. The experimental data are from [20]. 



TTT TTT 
Ji Jf 


A B[EX)expt 


B(i^A)calc 


4+^2+ 
6+^4+ 

9r ^ 6+ 
9r ^ 6+ 


2 4.3 ±0.3 

2 2.05 ± 0.03 

3 3.8 ±0.2 
3 8.0 ±1.3 


4.2 
1.9 
1.0 
8.2 


TABLE IL Calculated and 


experimental B{E2) values (in fm'*) for ^"^^Sb. The experimental data are from [28]. 


TTT JTV 


B(iJ2)cxpt 


S(i52)calc 


3+^4+ 
2+^3+ 
2~ 4+ 

22 > 32 

2+ 2+ 
2+ 3+ 

^2 "-"l 

2+ 4+ 

^2 ^ ^1 


36±11 
76±76 
< 26 

1 1 + 15.4 

37±30 
36±23 
8.4±4.5 


42 

53 

5.8 

0.85 

4.4 

5.4 

9.1 




FIG. 1. Experimental and calculated spectrum of ^*Cd. 






FIG. 2. Experimental and calculated spectrum of ^'^Ag. 





FIG. 3. Experimental and calculated spectrum of ^®Pd. 



FIG. 4. Experimental and calculated spectrum of Te. 



FIG. 5. Experimental and calculated spectrum of I. 



FIG. 6. Experimental and calculated spectrum of Sb. 
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